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Various  extensions  of  the  transmission  line  model  are  introduced  to  find 
the  resistance  for  current  flow  in  MOSFET  source/drain  regions.  The  ge¬ 
ometry  is  taken  to  be  a  rectangvilar  box  with  a  rectangular  contact  on  the 
upper  surface.  Explicit  formula  are  derived  by  aissuming  that  the  current 
flow  is  restricted  to  various  geometrical  planes.  Comparison  of  basic  re¬ 
sults  with  simulation  and  experimental  data  is  good.  Comparison  with 
simulation  results  for  misalignment  is  less  good. 
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Varioiis  models  have  been  suggested  in  attempts  to  find  resistance  formulae  for  cur¬ 
rent  flow  in  source/drain  regions  of  metal-oxide-semi  conductor  field-effect  transistors 
(MOSFET’s).  This  parasitic  resistance  becomes  relatively  more  important  under  de¬ 
vice  miniaturization,  and  its  dependence  on  geometrical  and  physical  parameters  is 
sought.  Measurement  of  resistance  across  single  contacts  ([l],[2])  indicate  large  vari¬ 
ability  neighbor  to  neighbor.  This  portends  the  use  of  statistics  in  estimating  the 
yield  of  a  device,  [3],  and  it  has  implications  for  the  quality  control  of  the  fabrica¬ 
tion  process.  Some  of  the  variability  is  remdom,  due,  say,  to  the  imperfect  removal 
of  photo-resist,  but  some  has  global  features.  Chief  of  the  latter  derive  from  mis¬ 
alignment:  contact  regions  have  been  etched  in  non-symmetrical  locations.  In  order 
to  estimate  the  effects  of  contact  size  and  location  some  two  and  three  dimensional 
models  are  constructed  here,  resulting  in  resistance  formulae  of  varying  complexity. 
Results  for  typical  parameter  ranges  are  compared  with  simulation  and  with  test 
data. 

The  simplest  model,  the  transmission  line  model  (TLM),  [4],  [5],  is  one-dimensional 
as  it  neglects  both  width  and  depth  variations.  Improvements  on  the  TLM  have 
been  provided  in  [6],  where  depth  dependence  was  included,  and  in  [7]- [9],  where 
modelling  of  contact  width  and  side  tabs  was  presented. 


The  physical  set-up  is  shown  in  Figure  1.  The  source/drain  region  is  taken  to  be  a 
rectangular  block  of  width  w,  depth  t.  The  contact  is  rectangular  of  length  b,  width 
Wc,  and  has  flange  distances  wi,W2,c  from  the  source/ drain  boundary,  as  shown. 

Current,  uniformly  distributed  across  the  end  surface  5,  flows  along  the  block  and 
exits  across  the  contact.  The  electrostatic  potential,  (f),  satisfles  a  zero  normal  deriva¬ 
tive  boundary  condition  on  bounding  surfaces  except  on  S,  and  except  on  the  contact 
where  an  ohmic  interfacial  resistance  holds.  This  gives  rise  to  the  boundary  condition 

(1.1)  <i>y  —  T(f>  =i0  on  the  contact 

where  r  =  lf{crpc). 

Here  cr  is  the  material  conductivity  and  pc  the  contact  resistivity. 

The  models  developed  in  this  paper  are  extensions  of  the  TLM,  and  of  the  work 
done  in  [6].  Model  0  aissumes  that  current  enters  the  contact  only  from  a  region 
directly  upstream  of  the  contact,  and  uses  the  TLM  formula.  Model  I  adds  to  this 
the  contribution  of  the  flanges.  Model  II  is  fully  three-dimensional,  extending  the 
analysis  in  [6].  The  key  assumption  is  to  prescribe  that  current  lines  remain  on 
planes  parallel  to  the  edges  of  the  box.  Current  flow  between  two  adjacent  planes 
is  then  two-dimensional,  and  the  full  result  is  achieved  by  integration  in  the  third 
direction.  Models  III  and  IV  are  variants  of  II  taking  into  account  different  current 
distributions. 
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Model  0  For  a  contact  of  width  Wc  the  series  or  front  resistance,  Rf,  is  given  by  the  TLM 

(2.1)  as  R,/Rf  =  (wc/L)ta.nh(b/L),  where  L  =  (tapc)^^^ 

and  Rt  =  pc/^^  is  the  sheet  resistance  of  the  material.  Substitution  of  sinh 
for  tanh  in  (2.1)  yeilds  the  end  resistance,  Rend-  Model  0  involves  neither  the 
flange  widths  wi,W2  nor  the  depth  t  (except  as  in  L).  The  models  introduced 
below  account  for  these  thickness  and  flange  effects.  As  a  result,  Model  0  cannot 
be  compared  properly  in  the  graphs  presented.  It  is  quoted  only  to  introduce 
the  formula  (2.1). 

Model  I  The  flanges  are  assumed  to  generate  additional  (parallel)  resistance,  computed 
by  viewing  the  flange  regions  as  having  width  t  and  depths  xu\,W2.  However, 
poor  results  from  this  at  small  values  of  Wc  (where  Rp  ^  is  expected)  made 
us  reflect  on  the  current  flow  interactions.  The  current  lines  squeezing  in  from 
the  flange  regions  are  competing  with  the  center  current  lines  which  have  depth 
t.  Detziils  of  the  current  flowing  in  from  the  sides  will  not  be  felt  by  the  contact 
until  that  ciurent  is  at  a  similar  distance  from  the  side  edge  of  the  contact. 
Consequently  (and  after  some  experimentation)  the  flange  regions  were  taken 
to  have  depths  where 


(2.2) 


<1  =  min(tyi,  t)  ,  t2  =  min(iy2,  t). 


The  front  resistance  is  then  given  by 


RJRf  =  iwcjL)  tanh(i/^)  +  (wi/Li)  tanh(u>c/Ti)  +  (W2/L2)  ta.nh{wc/ L2) 

(2.3)  where  Li  =  (t^apcY^^  ,  i  =  l,2. 

Model  II  This  is  an  attempt  to  construct  a  three-dimensional  approjcimation  to  current 
flow  in  the  complicated  geometry  of  Figure  1.  The  main  assumption  is  that 
current  lines  remain  on  planes  parallel  to  the  x-axis  (the  planes  are  not  par¬ 
allel  to  each  other).  The  planes,  and  the  reduction  of  the  three-dimensional 
problem  to  a  sequence  of  two-dimensional  ones,  are  found  by  the  following  two 
prescriptions. 

Prescription  (1):  Consider  the  end  surface  5  and  a  line  through  z  =  r.  The  line  divides  5 
into  two  areas.  See  Figure  2.  Consider  the  plaine  defined  by  extending  this 
line  parallel  to  the  x-axis.  The  assumption  that  current  lines  do  not  cross 
such  plames  implies  that  current  entering  area  A  exits  on  the  portion  of  the 
contact  on  lui  <  z  <  r.  The  geometrical  relation 

(2.4)  {r  —  wi)/wc  =  A/a,Tea.  o{  S  =  A/tw 

is  taken  to  define  the  plane  through  z  =  r.  It  follows  that  there  is  a  point  tq 
where  the  area  A  becomes  rectangular,  of  area  A  =  rot.  Hence  (2.4)  gives 

(2.5)  vq  =  wwi/(w  —  vuc). 
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It  is  now  permissible  to  concentrate  only  on  the  portion  of  the  source/drain 
region  defined  by  0  <  r  <  tq,  as  the  remaining  portion  will  have  the  same 
characteristics  under  this  model,  and  its  resistance  can  be  found  by  appro¬ 
priate  substitutions.  For  the  rectangle  defined  by  0  <  r  <  tq  the  plane 
passing  through  the  corner  y  =  t,  2  =  0,  occurs  at  r  =  r*,  say,  where 

(2.6)  r*  =  ‘lww\l{2'w  —  Wc). 

(See  Figure  3).  Planes  through  z  =  r  iot  r*  <  r  <  tq  intersect  the  bottom 
of  the  region,  giving  an  intercept  h{r)  and  a  plane  depth  ei(r)  where 


(2.7) 

h  =  [r(2u;  —  Wg)  —  2wwi]/we, 

(2.8) 

el  =  t^  +  4[tx;u;i  -  r(w  -  Wc)]^/wl, 

whilst  planes  through  2  =  r  for  wi  <  r  <  r*  intersect  the  side  of  the  region 
with  intercept  p(r)  and  plane  depth  e2(r)  given  by 

(2.9)  p^2t‘w{r  —  xv\)l{rwc), 

(2.10)  ej  =  r* -I- 4t*u;^[l  — 

Prescription  (2):  Consider  an  area  element  formed  between  the  planes  2  =  r  and  z  =  r dr. 

These  are  shown  in  Figure  4  for  the  cases  r*  <  r  <  tq  and  Wi  <  r  <  r*. 
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The  area  elements  are  replaced  by  the  sector  area  element  shown  in  Figure 
4(c),  of  equal  area.  Note  that,  from  (2.4), 


(2.11)  dA  —  {wt/ Wc)dr. 

The  current  flow  in  an  area  element  is  now  assumed  to  take  place  in  the 
sector  element  of  Figure  4(c),  and  the  equation  to  be  solved  is  Laplace’s 
equation  in  cylindrical  coordinates,  assuming  angular  independence.  The 
analysis  for  this  is  given  in  the  Appendix.  The  resistance  for  the  element 
dA  given  in  (.\.37)  is  integrated  over  S  to  give  Rp  as 


where  €1,62  are  given  in  (2.8).  Note  that  the  result  (2.12)  is  for  resistance 
in  the  region  0  <  r  <  tq.  If  the  contact  is  symmetrically  placed,  tq  =  w/2, 
and  R,/Rf  for  the  whole  device  will  be  twice  that  of  (2.12);  if  not,  the 
expression  (2.12)  has  to  include  a  similar  expression  obtained  by  replacing 
wi  by  u;2- 

Prescription  (1)  of  Model  II  implies  that  an  element  dr  of  the  contact  width 
is  proportional  to  dA,  an  element  of  5,  see  (2.11).  Since  current  is  uniformly 
distributed  across  S  this  means  that  equal  elements  of  the  contact  receive  equal 
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amounts  of  current;  that  is,  the  currrent  density  is  uniform  across  the  contact 


I 

1 

I 


width.  However  a  non-uniform  distribution  is  expected  as  current  entering  from 
the  flanges  will  be  confined  to  a  layer  close  to  the  sides.  .A.n  estimate  of  the  non¬ 
uniformity  is  provided  by  the  TLM.  that  is  by  Model  I.  The  current  density 
for  the  region  defined  the  contact,  i.e.,  ivi  <  r  <  tq,  is  taken  as  uniform.  The 
current  coming  from  the  flange,  0  <  r  <  u^i,  is  estimated  by  the  TLM.  as 
proportional  to 


(2,13) 

Combining  these  yields  the  replacement 


(2.14) 


dA  =  t 


1  -j-  ^cosh 


dr 


for  (2.11),  and  an  appropriate  change  in  the  formula  for  Rf  in  (2.12).  Also  the 
calculation  for  r*  is  changed  to  the  solution  to 


(2.15)  r-  =  2u;iA'(r*) 

where  K{r)  =  sinh(  °  ^){sinh 

L\  L\ 

Similar  changes  are  made  for  the  current  contribution  from  the  flange  u>2. 

Model  IV  The  introduction  in  Model  III  of  a  relation  between  dA  and  dr  which  is  non- 
uniform  also  giv'es  .4  as  a  function  of  r  (by  integrating  (2.14)  )  different  from 
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(2.4).  There  follows  new  formulae  for  £1,62 

=  +  Aw\K'^{r)  ,  e\  =  r-  +  4i^{l  -  — A'(r)}^ 

r 

which  are  incorporated  in  Model  IV. 

III.  Results  Various  results  derived  from  the  formula  (2.12)  and  the  corresponding  ones 
from  Models  III  and  IV  are  shown  in  Figures  7  to  10.  In  the  TLM  the  thickness  t 
is  parametrically  involved  only  in  the  definition  of  L  and  this  carries  over  to  most 
simulations  as  the  latter  are  derived  from  the  numerical  solution  of  two-dimensional 
partial  differential  equations  (in  r  =  0).  Loh  et  al,  [9],  have  carried  out  some  three- 
dimensional  simulations  and  show  differences  between  1-D,  2-D,  and  3-D  results 
for  jRfc,  the  Kelvin  resistance,  which  is  not  appropriate  for  the  model  constructed 
here.  There  are  more  simulations  for  the  2-D  case  and  we  use  these  for  comparison. 
However  we  note  that 

(1)  we  are  comparing  different  models,  and 

(2)  as  is  evident  from  formula  (2.12),  and  as  noted  in  [9],  Appendix,  the  thickness 
t  enters  the  3-D  models  as  a  parameter  to  be  chosen  in  relation  to  the  other 
length  scales  wi,Wc,L.  As  a  consequence  we  can  vary  t  to  see  its  effect  on  the 
results. 
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Figures  7  (a)  -  (d)  show  log(i?F/i?,)  versus  log(u;c/u;i)  for  L/iui  =  .358  for 
Models  I-IV  for  i/u^i  =  .1,  .5, 1, 10  respectively.  Simulation  data  from  [9],  Figure 
12,  is  also  shown.  It  is  noted  that  t/wi  =  .5,1,10  give  good  comparisons, 
indicating  some  insensitivity  to  t  over  this  range.  As  above  we  note  that  we  are 
comparing  different  models.  Figures  8  (a)  -  (d)  repeat  the  same  results  for  the 
case  Llw\  =  .567. 

Figure  9  shows  \og{REl Rs)  versus  log(u;c/u;i)  plus  simulations  from  [9],  Figure 
11,  for  f/ii/’i  =  1  and  for  Ljwx  =  .567  (upper  curve),  L/wx  =  .358  (lower  curve). 
These  show  reasonable  comparisons. 

Results  relevant  to  misalignment  estimates  are  shown  in  Figure  10,  where  also 
are  drawn  the  simulation  data  of  [10]. 

IV.  CQmmgntis 

There  are  small  differences  between  the  results  presented  here  and  the  simula¬ 
tions  in  Figures  7,8,9  and  substantial  ones  in  those  of  Figure  10.  Apart  from 
the  fact  that  Rend  is  more  volatile  than  Re,  there  are  significant  differences  in 
the  approaches,  which  bear  comment. 
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The  2-D  simulations  are  based  on  solving 


(4.1) 


=  <f>  on  the  contact, 


V2^)  =  0  off  it 

by  numerical  methods  over  the  plane  2  =  0,  imposing  a  constant  voltage  up¬ 
stream  at  a  station  such  as  S.  This  turns  out  to  be  equivalent  to  the  TLM. 
A  derivation  of  (4.1)  based  on  averages  across  a  thin  layer  is  given  in  [9].  .A.n 
alternative  is  as  follows.  Note  that  the  boundary-value  problem  posed  by  (1.1) 
is  equivalent  to  minimizing  the  integral 


(4.2)  7  =  /  +  r  I  6^dS. 

J  volume  J  contact 

Now  assume  at  small  t  that  the  volume  integral  may  be  approximated  by  t  times 
the  surface  integral  and  assume  that  the  r  term  applies  only  on  the  contact.  The 
zero  variation  of  J  yields  (4.1).  We  remark  in  passing  that  numerical  simulations 
based  on  (4.2)  may  be  more  efficient  than  3-D  finite  difference  schemes. 

The  work  in  [6]  improved  the  TLM  to  account  for  depth  effect  (to  first  order). 
It  showed 

(a)  that  (4.1)  must  be  supplemented  by  a  transition  layer  of  width  L  on  either 
side  of  the  contact  boundary.  A  more  complicated  boundary  value  problem 
holds  in  this  layer. 

(b)  that  the  effect  of  including  depth  is  substantial  on  RenDi  less  so  on  Rf. 
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The  formulae  presented  in  Section  II  do  not  include  the  depth  effects  calculated 
in  [6].  The  latter  could  be  included:  easily  in  Model  I;  approximately,  and  at 
the  expense  only  of  more  complicated  formula,  in  Models  II-IV.  However  the 
models  presented  do  take  into  account  the  fact  that  current  lines  have  a  three- 
dimensional  trajectory,  and  it  is  argued  that  they  may  be  more  realistic  than 
simulations  based  on  (4.1)  which  is  a  two-dimensional  approximation. 

The  divergences  between  the  two  approaches  as  seen  in  Figure  10  suggest  that 
the  modeling  is  inappropriate  for  misalignment  effects  in  one  or  both  the  ap¬ 
proaches.  The  modelling  here  assumes  current  flow  along  planes  parallel  to  the 
x-axis.  As  a  consequence  the  formulae  derived  involve  the  contact  length  b  as 
an  argument.  In  practice,  current  lines,  on  approach  to  a  square  contact,  may 
react  to  the  proximity  of  contact  side  edges  as  well  as  the  rear  edge,  thereby 
reducing  b  in  the  argument.  This  effect  shows  more  dramatically  in  Re.vd  which 
evaluates  sinh  of  arguments  than  in  Rf  which  evaluates  tanh. 
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APPENDIX 


We  wish  to  solve  Laplace’s  equation  in  cylindrical  coordinates,  with  no  angular  varia¬ 
tion.  That  is 

(.4.1)  d^<i>ldr^  -I-  (l/r)(9d/5r  -f-  d^<i>/dx'^  =  0 

for  the  boundary  value  problem  shown  in  Figure  5.  The  sector  occupies  the  region 
Ti  <  r  <  Tj  and  — oo  <  x  <  oo.  (The  source/drain  boundary  at  i  =  —(6  -f  c)  is  removed  in 
this  analysis  as  its  effect  on  the  results  is  assumed  to  be  weak.  An  approximation  which 
corrects  for  this  omission  can  be  introduced  without  difficulty).  The  sector  depth 

(.4.2)  e  =  r2  -  rj 

is  assumed  to  be  small,  and  an  asymptotic  result  based  on  e  <  1  similar  to  the  analysis  for 
the  plane  case  presented  in  [6]  is  obtained.  That  analysis  was  predicated  on  the  existence 
of  boundary  layer  solutions  in  regions  close  to  the  ends  of  the  contact,  i.e.  on  planes  i  =  0 
and  I  =  —b.  Consequently  there  are  five  regions  where  separate  solutions  are  found  (see 
Figure  6),  and  these  solutions  are  joined  in  the  conventional  approach  of  the  method  of 
matched  asymptotic  expansions  ([11]  or  [12]).  Solutions  for  the  five  regions  are  as  follows: 
Region  I:  This  region  contains  no  current  lines  and  the  potential  is  constant.  Hence 

(A.3)  4>  =  a  +  €T/3-^ - , 

where  a,  3  are  constants  to  be  found  subsequently. 
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Region.  II:  In  terms  of  scaled  variables  for  the  region  near  the  contact  end 
(*4.4)  m  =  (x  +  6)/e  ,  n  =  (r  -  ri)/e, 

an  asymptotic  expansion  for  the  solution  is 


(.4.5) 


^  =  a  +  €t[/3  +  n)]  + 


The  perturbation  potential,  (i)i,  satisfies  boundary  conditions 

on  n  =  0  ,  m  <  0 


(.4.6) 


d6ildn  =  0 
and  d<i>\/dn  =  1 


on  n  =  1  ,  all  m 
on  n  =  0,  m  >  0. 


In  the  plane  case,  [6],  an  e.xplicit  integral  e.xpression  was  found  for  4>i\  here  a  series  expan¬ 
sion  is  obtained. 

In  m  <  0 

(.4.7)  4>x  =  where  s  =  n  +  rje, 

N 

(A.S)  and  where  +  {l/s)F'^  +  N^F,\  =  0. 

Hence 

(.4.9)  F.v  =  .4,vJo(iVs)  -I-  B^YoiNs). 

The  boundary  conditions  that  Fv  =  0  on  n  =  0  and  n  =  1  determine  the  eigenvalues  N 
and  the  ratios  A/v  :  The  purposes  for  which  we  require  the  solution  do  not  require  us 
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to  complete  this  part  of  the  calculation. 

In  m  >  0 

The  condition  that  the  boundary  condition  (A. 6)  yields  the  correct  current  exiting  the 
contact  implies  that  <j>i  is  O(m^)  as  m  — *•  oo,  and  that 

(-4.10)  (l>i  =  -  s^/2]  +  k2lns  4-^:3  +  ^  e~^'^”'Gyj{s), 

iV 

(.4.11)  where  h  =  eri/{rl  -  rj),  ^2  =  Txrllt{r\  -  rj). 

The  solution  for  the  eigenfunctions  Gy  is  as  in  (-4.9)  with  constants  again  determined  (in 
ratio)  by  the  boundary  conditions.  The  continuity  of  61  and  (pim  across  m  =  0  is  required 
to  supply  sufficient  conditions  to  fix  the  constants  in  (.4.9)  and  similar  ones  in  (-4.10). 
There  remains  the  constant  ^3  in  (A.  10)  to  be  determined.  In  the  pleme  case,  [6],  this 
constant  was  determined  as  a  consequence  of  the  availability  of  an  exact  result.  However 
it  weis  noticed  there  that  the  constant  could  also  be  determined  from  the  physical  principal 
that  the  average  potential  increase  in  the  x-direction  is  a  consequence  of  current  flow  across 
the  contact.  This  latter  condition  is  invoked  here  and  gives 

(-4.12)  *3  =  (ri/4€A^)  {(1  +  ri/e)"*  -  (rj/e)^} 

-2(1  -f  rxitf  [2(1  +  ri/e)^/n(l  +  rj/e)  -  2(ri/€)^ln(ri/e)  -  A]} 

(-4.13)  where  A  =  (1  +  rj/e)^  —  (n/e)^ 
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Note  that  the  series  terms  in  (A.  10)  are  exponentially  decreasing  as  x  increases  into  the 
contact  region,  and  the  non-series  terms  carry  the  dominant  information  concerning  the 
current  flow  to  be  matched  with  the  solution  in  Region  III. 

Region  III:  The  scaling  relevant  here  is 

(.4.14)  ^  =  (r/e)^/2^,  r;  =  (r-ri)/e. 

The  analysis  is  a  little  more  sophisticated  than  that  for  Region  II,  so  we  present  more  of 
the  details.  The  differential  equation  (A.l)  becomes 

(.4. 15)  +  «^n/(7  +  ^/«)  =  , 

and  the  boundary  condition  (2.1)  becomes 

(.4.16) 

The  asymptotic  series  is  taken  ais 

(.4.17)  <^  =  <^o  d"  cr<^i  -f  ^t^4>2  -!-••• 

so  that  the  first  two  terms  satisfy 
(A. 18)  <^07iTj  +  +  n/e)  =  6 

with  boundary  conditions 

(.4.19)  <?!>On  =  0  at  rj  =  0  and  77  =  1, 
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(.4.20) 


and  dim  +  ^ii/Cl  +  '’■/«)  =  -dots  . 


with  boundary  conditions 

(A.21)  4>ir,  =  (f>Q  at  r;  =  0,  and  =  0  at  r;  =  1. 

The  solution  of  (A. 17)  and  (A. 18)  is 

(A.22)  <f>o  =  M{0 

and  that  of  (A.  19) 

(-4.23)  (f>i  =  M"{ri  +  ri/€)2/4  +  +  '*i/e)  +  -^(0. 

with  boundary  conditions  (A.21)  giving 

(A.24)  Ne/ri  -  iVrri/2e  =  M 

and  A’/(l  +  ^i/t)  ~  M"(l  +  rx/€)l2  =  0 

Hence,  eliminating  N, 

(A.25)  M"A  =  M2ri/e 

The  solution  of  (A. 25)  which  matches  with  the  leading  term  from  (A. 10)  for  Region  II 
gives 

(A. 26)  <t>()  =  a  cosh[/i(^  -I-  F)] 
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(A.27) 


where  F  =  ^  =  (2ri/eA)*^^. 


The  solution  for  <f)i  is  found  in  similar  fashion; 

(A.28)  <f>i  =  (1/2)q!^*{(1  +  rieyin{T]  +  ri/e)  -  (77  +  rxltf/2}  cosh7i(^  +  F) 

+k40i((  +  F)sinh/Lt(^  +  F)  +  (/?  +  a k^)  cosh  +  F), 

where 

Aki  =  +  ri/e)*  -  (ri/e)‘‘]/16  -  n{ri/e)^/4 

(■^•29)  -M^(l  +  ^'i/c)^[(l  +  rifeyin{l  +ri/e)  -  ri/eyin{ri/e)  -  A/2]4 

+(l/^^)(n/«)^(l  +  ri/£yin{ri/e) 

Refflon  IV:  The  scaling  used  near  the  contact  leading  edge  is 

(A.30)  p  =  x/e  ,  g  =  (r  -  ri)/e. 

Matching  with  the  solution  from  Region  III  indicates  that  a  series  in  powers  of 
required,  euid  the  solution  is  developed  as 

(A. 31)  <f>  =  acoshpF  +  {rey^^afip  sinh  pF  +  rec^jlP)?)  +  •  •  • 


Here  <^3  is  harmonic  and  satisfies  the  boundary  conditions 


(A.32) 


f 


d<f)2/dq  —  0  < 


on  g  =  0,  p  >  0 
on  g  =  1,  all  p, 


and  d(f>2/ dq  =  a  cosh  fiT  ong  =  0,  p<0. 


The  solution  for  (^2  is 


(.4.33)  <i>2  =  ak^T  sinh^r  +  {/3  +  oik^  +  afi'^[p^/2  —  {q  +  ri/e)^/^ 

+(1/2)(1  +  ri/e)^/n(9  +  ri/e)]}  cosh/iT  +  +  ri/e)  for  p  <  0, 

,v 

and 

<f>2  =  +  ok^)  coshp7  +  k^aV  sinhpF  4-  a  ^  +  ri/e)  for  p  >  0. 

N 

The  terms  Hs,  Tv  in.  the  infinite  series  above  are  linear  combinations  of  Bessel  functions 
and  are  evaluated  in  a  similar  fashion  to  those  in  the  Region  II  solutions.  However,  their 
exact  nature  is  not  required. 

Region  V:  The  asymptotic  behaviour  of  the  Region  IV  solution  as  p  — ^  00  provides  the 
solution  here.  This  is  found  to  be 

(.4.34)  6  =  (r/€)'^^iar  sinhpT  +  acoshpr+ 

r€[(^  +  04-3)  cosh pF  +  (l/2)fc4oF  sinhpF]  -}-••• 

The  first  term  in  (.4.34)  represents  the  uniformly  distributed  current. 

This  completes  the  analysis  for  generating  and  matching  the  solutions  in  Regions  I- 
V.  The  constants  a,0  eissumed  for  the  solution  in  Region  I,  equation  (A. 3),  may  now  be 
specified.  If  the  total  current  in  Region  V  is  I,  and  the  cross-sectional  area  is  S,  the  current 
density  is  I/S.  Equating  this  with  ad<l>/dx  from  (A. 34)  yields 

(.4.35)  a  =  (//5){,7(r/€)‘/VsinhpF}-T 
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(.-1.36) 


j  =  -%^(i-hAircoth/ir) 

■iri 


The  resistance  of  an  element  dS  may  now  be  calculated.  The  end  resistance  is  related 
to  the  potential  at  i  =  —6,  and  the  front  resistance  to  that  at  j  =  0.  Hence 

Re.\'D  ={ot  +  eTi3-\ - )/{jdS) 

(.4.37) 

=  {<T(r/e)^/^/isinh/irdS}“^ 

to  first  order.  Likewise  Rf  has  a  similar  form,  replacing  sinh  by  tanh. 
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Figures  7  (a)  -  (d~)  Graphs  of  log(i?f /i2j)  versus  log(u;c/u;i)  for  Models  I,  II,  III,  IV, 
for  t/wx  =  .1,, 5, 1,10,  respectively,  plus  simulation  data  from  [9].  The  graphs  may 
be  interpreted  as  keeping  t  and  loi  fixed  and  varying  Wc- 

Figures  S(a)  -  (d)  Same  as  Figure  7  with  L/wx  =  .567. 

Figure  9  Graph  of  log(i?£.vo/.Ri)  versus  log(u;c/u;i)  for  f/u-’i  =  1  and  for  Ljwx  =  .567 
(upper  curve),  L/wx  =  .358  (lower  curve). 

Figure  10  Graph  of  Re  versus  loi  to  show  the  effect  of  misalignment.  The  results 
are  for  a  square  contact  with  Wc  —  5/um  and  R,  =  20Q/D  and  for  t  =  .5,  l/zm.  The 
simulation  data  is  taken  from  [10],  Figure  4. 


Model  I 
Model  II 
Model  III 
Model  17 


O  O  O  O  Sloulatioa  Data 


L/W1*0 . 35a 


t/Wl*0 . 1 


FIGURE  7a 


0 

Log (fc/J 1 ) 


— 

yodel 

I 

yodel 

II 

______ 

yodel 

III 

yodel 

IV 

o  o 

Simula 

tion  Data 

L/Wl 

=  0 . 35S 

t:/Wl  =  0 

FIGURE  7b 


_  Model  I 

_  Model  II 

_  Model  III 

.  Model  IV 

O  O  O  O  Slaulatioa  Data 


L/W1=0.35S 

FIGURE  7c 


-1  0  1 

Log(Tc/T 1 ) 


-  Model  I 

-  Model  II 

-  Model  III 

.  Model  IV 

O  O  O  O  Sinulation  Data 


L/W1  =  0.35S  t:/Wl  =  10.O 


FIGURE  7d 


0 

Log(Tc/Tl) 


-  Model  I 

-  Model  II 

-  Model  III 

.  Model  17 

O  O  O  Siaulatioa  Data 


L/W1*0 . 567 


t/Wl 


FIGURE  8 


Log(Rf/Rs) 


O  O  O  O 


Model  I 


_  Model  II 

_  Model  III 


Model  IV 
Simulation  Data 


L/W1=0 . 567 


t/Wl=0 . 5 


FIGURE  Sb 


0 

Lag(Wc/W  1 ) 


1 


_  Model  I 

- ^ -  Model  II 

-  Model  III 

.  Model  17 

O  O  O  O  Siaulatioa  Data 


L/W1*0 . 567  .  O 


FIGURE  Sc 


0 

Log(Tc/T:) 


1 


o  o  o  o 


L/Wls 


2<Qdel  I 
lodel  II 
Jfodel  III 
Model  17 

Siaulation  Data 


.  567 


t/Wl-lO . O 


FIGURE  Sd 


FIGURE  9 


Model 

I  I 

wi  th 

ft 

II 

H 

• 

0 

Model 

I  I 

wi  th 

t  =  0  -  5 

Simulation  Data  -firom  C103 


FIGURE  lO 


